Prom reversing to hemispherical dynamos 
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We show that hemispherical dynamos can result from weak equatorial symmetry breaking of the 
flow in the interior of planets and stars. Using a model of spherical dynamo, we observe that the 
interaction between a dipolar and a quadrupolar mode can localize the magnetic field in only one 
hemisphere when the equatorial symmetry is broken. This process is shown to be related to the one 
that is responsible for reversals of the magnetic field. These seemingly very different behaviors are 
thus understood in a unified framework. 

PACS numbers: 47.65.+a, 52.65.Kj, 91.25.Cw 



Planets and stars seem to be able to generate hemi- 
spherical dynamos. During the Maunder minimum, the 
number of spots at the Sun's surface decreased and the 
sunspots were mostly localized in the southern hemi- 
sphere []J . It is believed that the magnetic field was then 
also dominant in this hemisphere [2|. Although Mars 
does not act as a dynamo anymore, it has been suggested 
that it did in the past and created an hemispherical mag- 
netic field Q. In addition to these observations, direct 
numerical simulations have clearly identified hemispher- 
ical dynamos Another surprising behavior is that of 
the Earth's magnetic field which reverses randomly in 
time. Reversals have been observed in various numer- 
ical simulations Q and in a laboratory experiment [f|. 
Among the numerical studies, events have also been re- 
ported where the magnetic field is much larger in one 
hemisphere Q. Reversing and hemispherically localized 
magnetic fields are two seemingly very different features 
of astrophysical, numerical and experimental dynamos. 

We explain these two different behaviors in a unified 
framework: they both result from the breaking of the 
equatorial symmetry of the flow. As direct numerical 
simulations have already exhibited the two behaviors at 
stake, we focus on a simple kinematic a 2 dynamo model 
in which equatorial symmetry can be broken. The ob- 
tained results are then generalized in the framework of 
dynamical systems theory: when the equatorial symme- 
try is weakly broken, reversals or hemispherical localiza- 
tion are shown to be two generic behaviors of any dy- 
namo which two most unstable modes are of dipolar and 
quadrupolar symmetries. 

We consider a sphere of radius R that contains a fluid 
of electrical conductivity a and magnetic permeability 
equals to that of vacuum /iq. The outside medium is 
vacuum. The inductive processes of the flow are modelled 
by a effects localized on two shells of vanishingly small 
radial thickness located at radius X-R an d £R- On the 
outer shell, an a effect converts the poloidal magnetic 
field into toroidal field. This latter is converted back into 
poloidal field on the inner shell. We thus write the a 
tensor in spherical coordinates (r, 0, </>) as 

otij = aeeS(r - S,R)8 ie 8 J e + a^8(r - xR)&i4>&j<t> ■ (1) 



The two scalars age and 
only and arc written 



are supposed to depend on 



= (Zi 
= a 2 



- bi cos(0) + ci cos(26>) + di cos(36») , 
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Odd dependences in latitude model an equatorially sym- 
metric flow whereas even ones break this symmetry. We 
solve the induction equation for this a 2 dynamo and look 
for axisymmetric eigenmodes B(r)e^' solutions of 



pB = Vx (orB) 



(3) 



for r < R and that match a potential field at r = R. 
Similar calculations were performed for a — uj dynamos 
with equatorially symmetric effects only [H, Q . From now 
on, we use the dimensionless growth rate p = p^aR 2 and 
a effects Ra\ = iaqg<x\, Rb± — fj,oabi... Moreover we set 
the position of the shells \ = 0-55 and £ = 0.9. 

We first consider an equatorially symmetric configu- 
ration with Rbi — Rb 2 — 20, Rd 2 = 25, all other pa- 
rameters being zero. Two modes are then linearly un- 
stable. The most unstable eigenmode is of dipolar sym- 
metry and is represented in fig. [1] (top). Its growth rate 
is p ~ 1.26. For this set of parameters, a second mode 
is unstable and its growth rate is only slightly smaller 
p ~ 1.21. This mode is represented in fig. [1] (bottom). 
It is even under equatorial reflection and is therefore of 
quadrupolar symmetry. The existence of two unstable 
modes depends on the values of the parameters. We ob- 
serve in general that the dipolar and quadrupolar modes 
have similar growth rates when the a effects are larger 
close to the poles. Roughly speaking, both hemispheres 
generate a magnetic structure. These structures are far 
from each other and thus only weakly coupled: the situ- 
ation with poloidal fields of same polarity, i.e. a dipole, 
does not strongly differ from a situation with opposite 
poloidal fields, i.e. a quadrupole. This argument, see 
for instance [1 01 ] , explains why a dipole and a quadrupole 
have similar growth rates and are likely to be unstable 
for close values of the parameters. 

Starting from this equatorially symmetric configura- 
tion, we now break equatorial symmetry by adding an 
a effect which has an even dependence on latitude. A 
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FIG. 1: Unstable magnetic field modes when the flow proper- 
ties are equatorially symmetric (see parameters in the text). 
The dashed circles are the two shells where the a effects are 
located. On the right hand side of the revolution axis (dot- 
ted), the lines of poloidal field are drawn. Lines of constant 
toroidal fields are shown on the left hand side. Top (resp. 
bottom): mode of dipolar (resp. quadrupolar) symmetry. 

first class of equatorial symmetry breaking leads to a 
decrease in the difference between the growth rates of 
the two modes, Ap. To illustrate this regime, we con- 
sider Rci = —Rc2 = P where P measures the symmetry 
breaking and Ra\ = Ra,2 = 0. When P increases, Ap 
decreases but both modes remain stationary until a bi- 
furcation occurs at P c = 0.78 where the eigenvalues of 
the two most unstable modes are equal. For P slightly 
larger than P c , they are complex conjugates. This corre- 
sponds to an exponentially growing oscillatory magnetic 
field. A snapshot of the time-evolution of the field during 
half a period is shown in fig. [5] for P = 1 . The dipole 
reverses and during a reversal, a quadrupolar structure 
is observed. The main effect of nonlinearities is to satu- 
rate the amplitude of the magnetic field (see discussion 
below): we have removed the exponential growth in this 
representation. We point out that the transition from 
stationary to oscillatory dynamo occurs for small symme- 
try breaking. For P = 1, the symmetry breaking effects 
(Rc\ and Rcz) are of the order of 5% of the equatorially 
symmetric terms. 

It has already been observed that oscillatory dynamos 



can be obtained from stationary modes when parameters 
are changed. Numerical studies of a — uj dynamos dis- 
play transitions from stationary to oscillatory decaying 
eigenmodes that most likely result from an interaction 
between two modes [111 ]. Competition between dipolar 
modes only has been studied for a 2 dynamos [l3 |: for 
particular forms of a effect, i.e. diagonal and dependent 
on r only, oscillatory modes can be observed. Here, we 
have shown that very weak equatorial symmetry breaking 
is sufficient to couple a dipolar and a quadrupolar mode 
so that a time periodic regime is obtained. When nonlin- 
ear terms are taken into account, this transition occurs 
through a saddle-node bifurcation and below the onset of 
bifurcation, turbulent fluctuations of the flow can trigger 
random reversals fl3j . Such time dependent fluctuations 
are not taken into account in our calculation. Never- 
theless, in the regime of random reversals, the snapshots 
obtained during a reversal would be similar to that dis- 
played in fig. O once a reversal has been initiated, its 
evolution is mainly deterministic and corresponds to half 
a period of the nonlinear oscillation obtained above the 
saddle-node bifurcation. 

We have also identified other terms that break equato- 
rial symmetry but do not lead to reversals: the two most 
unstable modes remain stationary and the difference be- 
tween their growth rates Ap increases. These modes dis- 
play a surprising spatial structure: the field is almost 
completely localized in only one hemisphere. This is il- 
lustrated with the set of parameters Ra\ = Ra2 = —P 
and Rci — Rc2 — and in fig. [3jthe most unstable mode 
is represented for P = 0.2. Note that the second most 
unstable mode is localized in the other hemisphere. To 
quantify this localization, we compute E n (resp. E s ) the 
energy of the field measured at the surface of the north- 
ern (resp. southern) hemisphere for a fixed radius r. We 
plot in fig. 2] the relative energy in the northern hemi- 
sphere fnh = En I (Pn + Es) as a function of P. For large 
P, the ratio f n h saturates to a limit value f°°. This 
value depends on r and tends to 1/2 when r becomes 
very large. Indeed, for large r/R, the field is dominated 
by the dipole for which f n h — 0.5. Nevertheless, we point 
out that even for r = 2R, 83% of the energy is still lo- 
calized in the southern hemisphere in the case shown in 
fig. [31 This effect can thus be observed at the surface of 
planets and stars even if dynamo action takes place much 
deeper in the fluid core. At first sight, one might believe 
that the localization is a straightforward consequence of 
the a-effect being more intense in one hemisphere. We 
present here the simplest form of symmetry breaking but 
we have obtained the same behavior with different ones, 
such as Ra\ = —Ra,2 = —P: the two a effects are en- 
hanced in different hemispheres but hemispherical local- 
ization still occurs. Even more surprisingly, we point out 
that a very weak symmetry breaking has a drastic effect 
on the field localization. In fig. [H for P > 0.05 i.e. as 
soon as Ra\/ Rb\ = 0.25%, more than 95% of the energy 
is located in the south hemisphere at r = R. 

The fact that a symmetry breaking of a few percents 



FIG. 2: Evolution of the magnetic field during a reversal obtained when the equatorial symmetry is weakly broken. From left 
to right and top to bottom, half a cycle is represented at time t = (arbitrary), t = T/10, 2T/10, 3T/10, 4.5T/10 and T/2, T 
being the period. The dipolar field evolves toward a quadrupole before a dipole of opposite polarity is recovered. 




FIG. 3: Most unstable magnetic mode for a symmetry break- 
ing that leads to field localization (see parameters in the text). 
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FIG. 4: Relative energy at the surface of the northern hemi- 
sphere f n h = E n /(E n + E a ) as a function of the intensity P 
of the equatorial symmetry breaking effect. The curves corre- 
spond to (thick) r — R, (dashed) r = 2R and (solid) r = 5R 
for the two unstable eigenmodes, the lower curve correspond- 
ing to the most unstable eigenmode. 



leads to such a strong modification of the geometry of 
the eigenmode is a consequence of the interaction be- 
tween two modes of instability. Indeed, hemispherical lo- 
calization can be understood as a competition between a 
dipolar mode, D(r) with amplitude d(t) and a quadrupo- 
lar mode Q(r) with amplitude q(t) from which the total 
magnetic field is determined, B = d(i)D(r) + g(t)Q(r). 
The amplitude equation at linear order for A = d + i q is 

A = (ti r + ini)A+(v r + iUi)A, (4) 

where the coefficients /i r , Hit v r and i>i are real. Searching 
for a phase equation, we write A — r exp (i0), and obtain 

= m + Vi cos (20) - v T sin (20) . (5) 

When the flow is equatorially symmetric, fii = v % = 
and u r measures the difference of growth rates between 
the two modes. When equatorial symmetry is broken, \Xi 
and Vi increase. If \Xi increases faster than i>i, a regime 
of reversals is obtained [l3| . 

Our results show that some terms that break equatorial 
symmetry do not lead to reversals, i.e. lead to Wj increas- 
ing faster than fii- This is responsible for hemispherical 
localization of the field. Indeed, if vi is large compared to 
Hi and v r , the solutions of ((5|) remain stationary and are 
~ 7r/4,37r/4 [tt]. These values correspond to a total 
field B oc D ± Q. Its spatial structure is easily under- 
stood: in one hemisphere the amplitudes of the dipolar 
and quadrupolar modes are of same sign and the field is 
strong, whereas in the other hemisphere, the dipole and 
quadrupole are roughly opposite and their sum results 
in a field of very low intensity. We have highlighted the 
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fact that, for hemispherical localization to occur, sym- 
metry breaking a effects do not need to be large. This 
traces back to the similar growth rates of the dipolar 
and quadrupolar modes. More precisely, the symmetry 
breaking effects, of order 2/3, need only to be larger than 
the growth rate difference v r . Both terms can thus re- 
main small. 

The role of nonlinear terms can be understood in this 
low dimensional approach. Terms in higher powers of 
A and A have to be taken into account in (fj|. Their 
main effect is to saturate the field amplitude r. They 
also modify the coefficients of ([5]), but the qualitative 
features of reversals are the same as those displayed in 
fig. [21 More precisely, nonlinear terms would change the 
critical value P c and the precise temporal dynamics of O. 
However, these two effects are of higher order in the field 
amplitude and remain small provided the modes are only 
slightly supercritical. This is also true for hemispherical 
dynamos. 

Thanks to the low dimensional approach we have 
shown that the results from the a 2 model are generic 
and can be used to describe the dynamics observed in 
astrophysical, numerical and experimental dynamos. Di- 
rect numerical simulations of convective dynamos have 
displayed hemispherically localized fields where "both 
quadrupolar and dipolar fields contribute nearly equal 
magnetic energy" even though "the convection flow has 
lost little of its equatorial symmetry in the presence of 
the hemispherical magnetic field" [4] . Numerical models 
of Mars have also obtained hemispherical dynamos [3(. 
Our results give a simple explanation for these localized 
fields: if the dipolar and quadrupolar modes have simi- 
lar growth rates, even a very weak equatorial symmetry 
breaking is enough to localize the field. As far as rever- 
sals are concerned, a numerical study of the Von Karman 
Sodium experiment displayed a transition from steady to 
oscillatory magnetic field resulting from the coupling be- 
tween a dipolar and a quadrupolar mode [l4l |. 

Studies aimed at modelling the Maunder minimum 
considered an oscillatory dipole in competition with an 
oscillatory quadrupole with similar pulsations [l5| : non- 
linear coupling terms between these modes can gener- 



ate hemispherically localized magnetic fields that remain 
time-periodic. The pure dipolar state is still a solution 
but is unstable towards a mixed polarity state. Further 
above onset, chaotic behaviors have also been observed 
flU IT(| . Our mechanism for localization also applies to 
this situation and does not require nonlinear effects. In 
contrast to the aforementionned scenario, it does not re- 
sult from a secondary instability of a pure dipolar mode 
but from the modification of the dipolar solution induced 
by the breaking of the flow's equatorial symmetry. 

We can speculate on the origin of the symmetry break- 
ing. First the flows are turbulent and velocity fluctua- 
tions are expected to break equatorial symmetry. Sec- 
ond, the boundary conditions can be asymmetric: the 
pattern of thermal flux at the core mantle boundary has 
been shown to influence the frequency of reversals [1 7] . 

Our simplified model captures some widespread fea- 
tures of dynamos. Since our results are constrained by 
symmetry considerations, more realistic models should 
display the same behaviors. In other words, all dynamos 
that have two unstable modes of dipolar and quadrupolar 
symmetry behave the following way. Breaking the flow's 
equatorial symmetry leads to two different regimes: a 
first class of symmetry breaking, associated to large /Xj 
in eq. @, results in reversing and oscillating dynamos. A 
second class, associated to large i/j, generates hemispher- 
ical dynamos. The fact that these apparently very differ- 
ent regimes result from the same flow feature is appeal- 
ing and strengthens the importance of the competition 
between dipolar and quadrupolar modes. In particular, 
we expect that several new observations of hemispherical 
magnetic fields will be reported in numerical, experimen- 
tal or astrophysical dynamos. 

Out of the dynamo context, our study shows that the 
competition between two modes of different symmetries 
leads to rich behaviors. It is known that reversals can 
be observed in large scale flows generated over a tur- 
bulent background in thermal convection or in periodi- 
cally driven flows [l8| . We expect that a weak symmetry 
breaking of the forcing can also generate a localization 
of these large scale flows. This could strongly affect the 
mixing properties or the thermal transfer efficiency of 
these systems. 
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